Defect Formation Preempts Dynamical Symmetry Breaking in Closed Quantum 

Systems 



Carmine Ortix, 1 Jorrit Rijnbeek, 2 and Jeroen van den Brink 1 

institute for Theoretical Solid State Physics, IFW-Dresden, D01171 Dresden, Germany 
2 Insitute-Lorentz for Theoretical Physics, Universiteit Leiden, 2300 RA Leiden, The Netherlands. 

(Dated: January 15, 2013) 

We show that no matter how slowly a quantum-to-classical symmetry breaking process is driven, 
the adiabatic limit can never be reached in a macroscopic body. Massive defect formation preempts 
an adiabatic quantum-classical crossover and triggers the appearance of a symmetric non-equilibrium 
state that recursively collapses into the classical state, breaking the symmetry at punctured times. 
The presence of this state allows the quantum-classical transition to be investigated and controlled in 
mesoscopic devices by supplying externally the proper dynamical symmetry breaking perturbation. 
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Introduction. - The relation between quantum physics 
at microscopic scales and the classical behavior of macro- 
scopic bodies has been debated since the inception of 
quantum theory. The fundamental difference is that 
while in quantum mechanics all configurations equivalent 
by symmetry have the same status, in classical physics 
one of them is singled out - the symmetry is broken. This 
spontaneous symmetry breaking causes a macroscopic 
body under equilibrium conditions to have less symme- 
try than its microscopic building blocks [1 . Supercon- 
ductors, antiferromagnets, liquid crystals, Bose-Einstein 
condensates and crystals all exhibit spontaneously bro- 
ken symmetries. The general idea is that when the num- 
ber of microscopic quantum constituents, which depend- 
ing on the system corresponds to the number of Cooper 
pairs, particles or spins, goes to infinity the quantum 
system undergoes a phase transition into a state that 
violates the microscopic symmetries. From a purely the- 
oretical perspective, spontaneous symmetry breaking is 
thus related to a singularity of the thermodynamic limit. 
If this thermodynamic singularity is present, a symme- 
try broken groundstate exists. This observation can be 
formalized into an exact statement on the existence of a 
broken symmetry state 013], but it makes no assertion on 
whether or how it can evolve out of the symmetric state, 
nor on the dynamics of quantum-to-classical transitions. 
So the question remains how a continuous symmetry can 
be broken dynamically. 

We investigate this by considering a symmetry break- 
ing field that slowly drives an arbitrary large but fi- 
nite system from quantum to classical. If the symme- 
try breaking process were fully adiabatic, the effect of 
the driving would correspond to subjecting the system 
to a quasi-static symmetry breaking field. In this case 
the timescale at which the classical state is singled out 
and symmetry is broken becomes shorter and shorter as 
the size of the system grows. However, that the time- 
evolution be adiabatic is not evident. The adiabatic the- 
orem states that under slow enough external perturba- 
tions, there are no transitions between different energy 



levels. When the distance between energy levels is ex- 
ponentially small, the adiabatic evolution is hampered 
and transitions between levels become unavoidable. Such 
transitions and their associated defects can in princi- 
ple strongly affect the symmetry breaking process. Here 
we show that no matter how slowly a symmetry break- 
ing perturbation is driven, the adiabatic limit cannot 
be reached. Defect formation turns out to be so per- 
vasive that it preempts an adiabatic quantum-classical 
crossover in macroscopic systems. The existence of this 
non-adiabatic regime is consistent with the recent dis- 
covery [4] that adiabatic processes in low dimensional 
systems with broken continuous symmetries are absent. 

The far-from-equilibrium time-evolution caused by a 
symmetry breaking field has remarkable consequences. 
We will show that in any large finite system the non- 
equilibrium state does not break the symmetry. However, 
it recursively collapses into the purely classical state: it 
breaks the symmetry at punctured times, resulting in a 
Dirac comb of symmetry broken, classical states. This 
Dirac comb of quantum-classical transitions can be in- 
vestigated in mesoscopic devices by supplying a proper 
dynamical symmetry breaking perturbation. 

Spontaneous symmetry breaking. - Even if a priori 
spontaneous symmetry breaking is an intractable prob- 
lem involving a near infinity of interacting quantum de- 
grees of freedom, there is a representative, integrable 
model that exhibits spontaneous symmetry breaking: the 
Lieb-Mattis model [5] ■ It is the effective collective Hamil- 
tonian that underlies the breaking of the SU(2) spin rota- 
tion symmetry in generic Heisenberg models with short- 
range interactions. Very similar collective models under- 
lie the breaking of other continuous symmetries as the 
gauge invariance in superconductors or the translational 
symmetry in quantum crystals 0. The results that we 
will present here are therefore robust and generic. 

The quantum-classical symmetry-breaking transition 
is manifest in the Lieb-Mattis model once a symmetry 
breaking field H, in this case a staggered magnetic field, 
is introduced. Before turning to dynamical symmetry 
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breaking we first summarize a few essential features of the 
Lieb-Mattis Hamiltonian. The Hamiltonian is defined for 
spins 1/2 on a bipartite lattice with sublattices A and B, 
where Sa,b is the total spin on the A/B sublattice with 
z projection S A ^ B : 

H lm = 2 ^Sa-Sb-H(S* a -S b ). (1) 

Every spin on sublattice A interacts with all spins on 
sublattice B and vice versa with an interaction strength 
2\J\/N (which depends upon the number of sites N). 
Taking H — 0, the model can be solved by introducing 
the total spin operator S = S,4 + S#. The eigenstates 
of the Hamiltonian are then \Sa, Sb, S, M) where S, M 
indicate the total spin and its z axis projection, whereas 
Sa,b are the total sublattice spin quantum numbers. The 
ground state is symmetric and corresponds to an over- 
all S = singlet with Sa, b maximally polarized and is 
characterized by zero staggered magnetization. Magnon 
excitations carry an energy J and are realized lowering 
either Sa or Sb- The S ^ quantum numbers label a 
tower of states with energy scale -Ethi n = J/N, which is 
also referred to as the thin spectrum. It is thin because it 
contains states that are so sparse and of such low energy 
that their contribution to thermodynamic quantities van- 
ish in the thermodynamic limit [6]. Nevertheless, when 
N — > oo the thin spectrum excitations collapse and form 
a degenerate continuum of states. Within this contin- 
uum, even an infinitesimally small symmetry breaking 
perturbation H is enough to stabilize the fully ordered 
symmetry broken ground state - the system is inferred 
to spontaneously break its symmetry. The finite symme- 
try breaking field H couples the thin spectrum states so 
that the eigenstates |n) = X)s u sl^) °f the Lieb-Mattis 
model become wave packets of total spin states. In the 
continuum limit, where N is large and <C S <C N , the 
corresponding low-energy effective Hamiltonian is [3 



H = 



HN 



n 2 



N ' 



(2) 



where n is the conjugate momentum of the total spin 
S. The eigenstates Ug are harmonic oscillator states of 
order n, with n odd in order to meet the boundary con- 
dition S > 0. One can easily show the singular nature of 
the thermodynamic limit in the n = 1 ground state by 
calculating the expectation value of the order parameter 
2 (S A — S B ) ~ N e~" s where the dimensionless parame- 
ter w 5 = N^ 1 ^/4J/H. When sending first H -> and 
then N — > oo, the singlet state appears as the ground 
state, which respects the spin rotational symmetry, i.e. 
2(S A — S B ) = 0. Taking the limits in opposite order, 
one finds that the ground state corresponds to the fully 
polarized antiferromagnetic Neel state with a fully de- 
veloped order parameter 2{S A — S Z B ) = N. In this case 
Ethin = V J H represents the typical energy of the excita- 
tions labeled by n that now act as a dual thin spectrum. 
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FIG. 1: The three different regimes for the behavior of the 
density of defects in the t — to plane. Times have been mea- 
sured in unit of the freeze out time t. The bold lines indicate 
the crossover among the different regimes whereas the straight 
line limits the physical region with t > to- The inset shows the 
setup of the symmetry breaking field considered throughout 
this work. 



Adiabatic-impulse approach. - Let us now consider the 
dynamical case and turn on the symmetry breaking field 
linearly in time H (t) = St, with ramp rate r5. At initial 
time to we start out with a field H(to) = H [see the 
inset of Fig. [I] and the wavefunction of the system cor- 
responding to this static ground state. We introduce Hq 
in order to have a cutoff that guarantees the continuity 
of the wavefunction basis. Lateron we will consider the 
limit H — > which corresponds to an initial state that is 
a completely symmetric singlet. To capture the dynam- 
ics of the symmetry-breaking transition we first use the 
quantum Kibble-Zurek (KZ) theory The essence 

of the KZ theory of non-equilibrium phase transitions [9 
is a splitting of the dynamics into a nearly critical im- 
pulse regime where the system's state is effectively frozen 
and a quasi-adiabatic regime far from the critical point 
where transitions among the instantaneous eigenstates 
of the Hamiltonian arc neglected. This splitting defines 
the so-called adiabatic-impulse approximation |10j . In 
particular, the critical impulse regime occurs whenever 
the characteristic relaxation time r(t) — h/ E t hin(t) is 
much larger than the timescale t on which the Hamilto- 
nian is changed. On the contrary for r(i) << t, the sys- 
tem's state is able to adjust to the changing symmetry 
breaking field, and the transitions among the dual thin 
spectrum excitations can be neglected. The crossover 
between the two regimes is determined by Zurek's equa- 
tion [7] r{t) = t and defines the freeze-out time 
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For an initial time t 3> t, the system's dynamics will 
thus be nearly adiabatic [c.f. Fig. pH. Strictly speaking, 
in the true adiabatic limit (to/t — > oo) the probability of 
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FIG. 2: (a) Asymptotic value of the density of defects as a 
function of the initial time to over freeze-out time t. (b) Time 
evolution of the density of defects for different values of the 
initial time to. Times are measured in units of the freeze-out 
time t. The curves are independent of N . 
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switching thin spectrum levels will be vanishingly small. 
To quantify this, we calculate the fidelity of the snapshot 
ground-state wavefunction [TT] f{t) — \ (ug(t) | ifj{t)) | 2 
with ip(t) the actual ground-state wavefunction and the 
associated density of defects T>(t) = 1 —f(t) [c.f. Fig. [2]. 
The ramp rate of a nearly defect-free quench is seen to be 
bounded by S <C \J Hq J / h. Therefore the limits 5 — >• 
and Hq — > do not commute. In other words: no matter 
how slowly one drives the symmetry breaking field, if the 
initial symmetry breaking field Hq is sufficiently small, 
the adiabatic limit can never be reached. 

Besides just an adiabatic time-evolution, the KZ anal- 
ysis renders two non-trivial regimes for the dynamics 
of an initially symmetric ground state [c.f. Fig. 
First the evolution takes place in the impulse regime 
(to <C t <C t ) where the spin singlet state is effectively 
frozen and changes only by a trivial overall phase factor. 
In this regime, the density of dynamically generated de- 
fects grows continuously in time [c.f. Fig. j2^b)] . This 
evolution lasts until t > t. At the freeze-out time t, de- 
fect formation stops and the defect density saturates [c.f. 

Fig. 2] to V(i) - 1-8 (t /t ) 3/4 . The defect density thus 
tends arbitrarily close to 1 for small enough initial sym- 
metry breaking field Ho- In this case the system actually 
reaches a state that is a superposition of an exponen- 
tially large number of thin spectrum excitations. The 
quantum-classical crossover is thus accompanied by mas- 
sive defect formation, which in the end preempts such 
a crossover. In the subsequent adiabatic regime t t 
no new excitations are created. The corresponding evo- 
lution of the wavefunction can be analyzed as follows: 
we first expand the frozen initial ground state as a su- 
perposition of the dual thin spectrum eigenstates at the 
freeze out time t as Ug(to) = J2 n c nU 7 g(t), where the co- 
efficients c n are non zero only for odd values of the quan- 
tum number n. Since the evolution is considered to be 
adiabatic, we may write the time evolution of the wave- 
function as ips(t) — En c « M s(0 e l7 "(*) e~ ia "( t ', where 
7„(£) is the Berry phase [H] associated to the eigenstate 
In) which turns out to be zero and we have defined the 



FIG. 3: (a) Time-dependence of the real part of the dimen- 
sionless parameter to for different values of the initial time to 
for N — 10 2 spins and a freeze-out time t = 1. All times 
are measured in units of 4 J/8, (b) Same for a fixed initial 
time to = 10~ 2 and different values of the number of sites iV. 
(c),(d) Same as panels (a) ,(b) for the imaginary part of the 
dimensionless parameter ui. 
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the time increases, the various thin spectrum eigenstates 
all pick up a different dynamical phase leading to quan- 
tum interference. However for t/. = (1 + ^|^) 2 / 3 i with 
k integer, the interference is fully constructive and the 
wavefunction corresponds to the instantaneous ground 
state u\ (t to 1 1) as easily follows by considering that the 
coefficients c„ depends on the ratio to/t alone. The 
system's state then corresponds precisely to the snap- 
shot ground state of a symmetry broken Lieb-Mattis 
model subject to a renormalized staggered magnetic field 
H R = H x H /[h 2 S 2 /{2J)} 1 / 3 . Thus when the initial 
symmetry breaking field Ho vanishes, the symmetric sin- 
glet state becomes fact at any recursion time. 

Exact quantum theory. - A full description of the in- 
terference effects in a highly non-adiabatic state within 
the adiabatic-impulse method is in practice impossi- 
ble. We can, however, explicitly monitor quantum 
phase interference effects by constructing the exact 
non-equilibrium wavefunction for Hamiltonian Eq. ||2j). 
Within the Feynman path integral approach, it can be 
shown |13| that the propagator has a spectral decompo- 
sition Q(S B , t B \S A , t A ) = £„ ^ n s *(t A ) [t B ) in terms 
of a complete set of wavefunctions of the form 
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where the quantal phase 4>(i) and the complex dimen- 
sionless parameter ui(t) are uniquely determined by the 
classical Euler-Lagrange equation of motion whereas the 
quantum number n takes only odd values in order to 
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FIG. 4: Asymptotic behavior of the real (a) and the imaginary 

(b) part of the dimensionless parameter ui in the limit to — > 0. 

(c) Asymptotic behavior of the order parameter as a function 
of time measured in units of the freeze-out time t. 



meet the boundary condition S > 0. Different sets of 
wavefunction of the form Eq. ^ exist and correspond 
to take different pairs of linearly independent solutions 
to the equation of motion. This allows us to choose two 
particular solutions guaranteeing that at the initial time, 
^lj(to) = u s(to) so that the exact wavefunction remains 
an n = 1 state of the form Eq. Q. The wavefunction 
is thus completely characterized by a dimensionless pa- 
rameter u> (which plays a similar role as cog in the static 
case) that now has both a real and an imaginary part. 
The resulting time dependence of the real part of lo for 
different values of to is shown in Fig. [3|a) . By decreas- 
ing the initial time, it develops a series of sharp peaks 
which eventually leads to a Dirac comb structure in the 

to — > limit [c.f. Fig. [4|a)] with singularities at instants 

a 

~ t(|fc7r+ ^f 5 -) 3 ■ The imaginary part of lo has a 
time dependence as in Fig. [3jc) . As the initial time de- 
creases, it approaches a characteristic tangent-like behav- 
ior [c.f. Fig. H (b)] with singularities appearing precisely 
at the Dirac deltas of Re(u>). The limiting behavior of 
lo is universal since it scales with iV -1 as the number of 
sites is varied [c.f. Fig. |3jb),(d)]. 

The exact time development reveals that when the 
symmetric system is subject to the symmetry breaking 
field, a non-equilibrium state forms that is intermediate 
between a pure quantum symmetric and a pure classical 
state. It is a vast superposition of S =/= thin spec- 
trum excitations, with complex amplitudes. This state 
does not break the SU(2) symmetry, as a computation 



of (S\ — S%) directly demonstrates. As time evolves, 
this non-equilibrium state develops smoothly, until at a 
certain moment the system's state corresponds precisely 
to the classical symmetry broken ground state - the fully 
polarized Neel ground state of the symmetry broken Lieb- 
Mattis model. This classical state forms at punctured 



times t 1 /. ~ t (|fc 7r + ^p) 3 . At any other instant, the 
spin rotation symmetry is restored. This is in agreement 
with the adiabatic-impulse analysis which does not al- 
low symmetry breaking of a symmetric state when quan- 
tum phase interference effects are neglected. As a result, 
the time evolution of the order parameter is character- 
ized by a comb structure 2 (S% — S%) = N J2k>o ^t.t 1 
[see Fig. |4jc)] which corresponds to the periodic emer- 
gence of the symmetry broken states at punctured times. 
These instants are related to the freeze-out time alone, 
indicating the non-equilibrium nature of this dynamical 
symmetry breaking phenomenon. The freeze-out time 
can be experimentally tuned by changing the ramp rate 
of the symmetry breaking field and a quantum-classical 
transition can be induced in individual mesoscopic quan- 
tum objects by supplying a proper dynamical symmetry 
breaking perturbation. In the case of an infinitely sud- 
den quench (5 — ¥ oo), the freeze-out time vanishes and 
the punctured times of symmetry broken classical states 
collapse onto each other. In the contrary, asymptotically 
adiabatic limit (S — > 0), the first punctured time of sym- 
metry broken state diverges: the system never breaks its 
symmetry. 

Conclusions. - In the dynamical realm the quantum- 
classical symmetry breaking transition is thus character- 
ized by far-from-equilibrium processes. The exact theory 
shows that no matter how slowly the symmetry break- 
ing process is driven, defect formation prevents an adi- 
abatic quantum-to-classical time evolution. In a closed 
system, therefore, a stable symmetry-broken state can- 
not evolve out of a symmetric quantum state - neither 
spontaneously nor by driving it. 
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